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We have investigated how helium atoms are distributed 
within a mixed 3 Hejv 3 - 4 Hejv 4 large drop with N3 S> N4. For 
drops doped with a SF6 molecule or a Xe atom, we have found 
that the number of 3 He atoms within the volume containing 
the first two solvation shells increases when N4 decreases in 
a way such that these dopants may be in a superfiuid en- 
vironment for N4 > 60, which gradually disappears as N4 
decreases. The result is in qualitative agreement with recent 
experimental data. 
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In a recent experiment, Grebenev et al have carried 
out the equivalent of the Andronikashvili experiment |^] 
in a microscopic system, namelly a mixed 3 He- He drop 
consisting of about 10 4 atoms doped with an oxigen car- 
bon sulfide (OCS) molecule. By analyzing the infrared 
spectrum of OCS, Grebenev et al (see also Ref. ||) con- 
clude that the molecule freely rotates when a number of 
He atoms large enough coats the impurity preventing 
the 3 He atoms, which are in the normal phase at a tem- 
perature of the order of 150 mK |QJj|, from getting too 
close to the OCS molecule. That number is of the orden 
of 60, in excelent agreement with path integral || and 
variational Monte Carlo calculations. It is remarkable 
that the presence of the impurity, which causes the He 
density to rise up to several times the saturation value, is 
not destroying its superfiuid character, and that in spite 
of the high densities reached, the first solvation shell re- 
mains liquid [|| . An indication of this fluidlike behavior is 
that the peak density in the first solvation shell continues 
to increase as the second shell grows ||. 

Even if the intrepretation of the microscopic An- 
dronikashvili experiment is on a firm basis, a remaining 
major question is how 3 He is distributed around the 4 He- 
plus-impurity complex, and in general, how liquid 3 He 
is dissolved into 4 He droplets at very low temperatures. 
These are the questions we want to address in this work. 

At zero temperature, it is known that the maximum 
solubility of 3 He in the bulk of 4 He is ~ 6.6% fuj. For 
liquid 4 He systems having a free surface, it is also known 
that a large amount of 3 He is accumulated on the free 
surface occupying Andreev states [jlllH before it starts 
being dissolved into the bulk. In the case of drops made 
of up to several thousands of atoms, the surface region 
constitutes a sizeable part of the system fl3|| , and the 
surface has a large capacity of storing 3 He atoms before 
they get inside the drop [[l4| . Due to the wide free surface 



of both isotopes Jl5|,|T6| and to the low surface tension of 
the 3 He- 4 He liquid interface fill , one expects that this 
region plays a prominent role when it constitutes a large 
part of the system or, as in the present case, when it is 
close to the foreing atom or molecule. 

The structure and energetics of mixed, doped or not, 
helium droplets has been addressed using a finite-range 
density functional [pd[ . That work was carried out be- 
fore the experiments reported in Ref. (IJ] , and the empha- 
sis was put on improving the density functional to bet- 
ter describe the thermodynamical properties of the liquid 
mixture, and to study rather small mixed droplets with 
N± 3> A3. Our main goal here is to apply the density 
functional method to droplets whose characteristics are 
closer to those of the experiments, with the restriction of 
spherical symmetry for the He-impurity potential for the 
sake of simplicity. We have considered Xe and SF6 as 
dopants, using for the later a spherically averaged inter- 
action potential. The Xe-He potential is weaker than the 
SF 6 -He one. In this respect, our results for that atomic 
impurity should better represent the experimental ones 
for OCS even if this linear molecule produces deforma- 
tions in the helium drop that we have not considered 
here. The density functional and the treatment of the 
impurity are thoroughly described in Ref. Jlif . 

The large number of 3 He atoms in the droplets we 
are describing (A3 > 1000) allows us to employ an Ex- 
tended Thomas-Fermi method to describe the fermionic 
component of the mixture. We have used for the 3 He 
kinetic energy density the expression given in Ref. [fl8|| , 
which contains up to second order density gradient cor- 
rections to the standard ~ p^ 3 (r) expression, where p$ 
and P4 will denote the particle density of each isotope. 
We have checked that this density functional reproduces 
accurately the Hartree-Fock results |lj] obtained for the 
largest drops there studied (see also Refs. ]5|Jl^|). 

Figure |l| displays the situation in which a 4 He728 drop, 
whose size is large enough to clearly distinguish in it a 
surface and a bulk region, is coated with an increasing 
number of 3 He atoms, and the limiting situation of the 
same drop immersed into liquid 3 He. The evolution with 
A3 of the 3 He concentration inside the 4 He drop, defined 
as X3 = P3/(p4 + p?,)\buik, is shown in the insert. Several 
interesting features of this figure are worth to comment. 
First of all, a fairly large amount of 3 He is needed before 
it is appreciably dissolved in the bulk: for A3 = 1000, 
P3 near the origin is^l.4xlO~ 8 A . The solubility is 
appreciably reduced by finite size effects. Indeed, one can 
see from the insert that the limiting solubility into the N4 
= 728 drop is ~ 2.5%, as compared to the 6.6 % value in 
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the liquid mixture. It is also worth to see that for large 
N 3 droplets, the bulk solubility is slightly higher than the 
limiting solubility, indicating that finite size effects still 
appear in rather large drops. Another manifestation of a 
finite size effect is that the average 3 He density is above 
the saturation value even for the larger drops, showing 
that the existence of the outer 3 He surface still causes a 
visible density compression. 

Due to the high incompressibility of helium, the bulk 
density of 4 He decreases when 3 He is dissolved, and the 
rms radius of the 4 He drop manifests a peculiar N 3 behav- 
ior. It decreases when N3 increases up to a few hundreds 
due to the initial compression of the outermost He sur- 
face, and then steady increases as 4 He is pushed off the 
center by intruder 3 He atoms. This is a very tiny effect 
anyway. For example, we have found that the rms radius 
of the 4 He728 drop is 15.70 A. It decreases when 3 He is 
added, reaching a minimum value of 15.64 A for N3 ~ 
250, and then it steadily increases up to 16.11 A for N 3 
= 10000. The rms radius of the 4 He728 drop immersed 
into liquid 3 He is 16.14 A. 

When a SF 6 molecule is captured by a helium drop, it 
moves into the bulk producing a drastic rearrangement 
of the drop density around it pfj|-p3[|. For large 4 He 
droplets, it is especially noteworthy the appearence of 
two high density solvation shells with a density-depleted 
region in between. It is then natural to ask about the pos- 
sible existence of Andreev-like states arising at the 'inner 
4 He surface', and whether a large number of 3 He atoms 
can be stored there, producing an onion-like structure of 
alternative 4 He and 3 He shells around the impurity, or 
even if when A3 ^> N4, the later can displace the former 
in the first solvation shell. 

Figure || shows the density profiles of several 4 He728 + 
3 HeAT 3 + SF 6 droplets, giving a positive answer to the 
first question, and a negative answer to the other two. 
We have found that indeed, about one 3 He atom is in 
the inner surface, but that 3 He is mostly coating the 
4 He-plus-impurity complex, as in undoped droplets. To 
check this result we have started the calculations from 
different initial shapes, some having the onion-like form 
mentioned before. It has turned out that these always are 
high energy, metastable configurations, and the mixed 
droplet eventually evolves towards stable configurations 
of the kind shown in Fig. ^. The larger zero point motion 
energy of 3 He makes it energetically more advantageous 
to fill the first solvation shell with 4 He atoms, and 3 He is 
expelled to the outer region of the drop. 

We are now in position to discuss a physical situation 
relevant to the microscopic Andronikashvili experiment. 
We first observe that the first solvation shell |l9| can host 
~ 23 4 He atoms in the case of SF$ as a dopant, and ~ 
15 atoms in the case of Xe ]20|-p2|]. According to Refs. 

these numbers are too small for the 4 He droplet be- 
ing superfluid. It is thus crucial to know how the second 
solvation shell is built, especially what is its composition. 
Too many 3 He atoms in that shell might shrink or even 
wash out the superfluid environment around the dopant. 



The density functional method cannot tell whether a 
given configuration is superfluid or not, but it can give a 
quantitative answer to its local composition because it is 
able to reproduce available microscopic density profiles 
pOlpTpi] . We present examples of such compositions in 
Figs. H and f|. 

Figure m shows the density profiles for 4 HeAr 4 + 3 He 1000 
+ SF 6 and 4 HeAr 4 + 3 Hei 00 o + Xe with N 4 = 35, 60 and 
100. We have carried out calculations for two different 
dopants to see what is the influence of the He-impurity 
potential on the results. It turns out that a weaker at- 
tractive potential favors the mixing of both isotopes in 
the whole allowed volume (the Xe-He and SFg-He poten- 
tials are plotted in Ref. p^] , for instance). However, this 
is in part a first glance effect, since the number of 3 He 
atoms in the first solvation shell around Xe is less than 
one (see Fig. [|). Rather, the relevance of Fig. || lies in 
that it shows how 3 He is filling the second solvation shell 
as N4 decreases. 

A more quantitative look at this phenomenon is pre- 
sented in Fig. ||, where we have plotted the number of 
atoms of each isotope as a function of the radial distance 
to the center of the drop. Notice that for a given im- 
purity, the number of 4 He atoms in the first solvation 
shell (extending up to ~ 5.5 A) is sensibly the same for 
the three selected N4 values. It is also worth to look at 
the ratios N 3 /(N4 + A3) within the second solvation shell 
which extends from ~ 5.5 to ~ 8.5 A. In the SFg case, 
they are ~ 8% for N 4 = 100, - 29% for A 4 = 60, and 
~ 65% for N4 — 35. Considering the content of the two 
shells, these ratios are ~ 5%, ~ 19%, and ~ 41% which 
correspond, respectively, to 3, 10, and 22 3 He atoms. The 
values for Xe are slightly smaller. These numbers make 
it quite plausible that a SFg molecule or a Xe atom in 
a 4 Hejv 4 + 3 Heiooo drop is in a superfluid environment 
when N4 = 100 or 60, whereas it is not when N4 = 35, 
as the microscopic Andronikashvili experiment indicates 
for OCS. 
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for useful discussions. This work has been performed un- 
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alitat of Catalunya. 
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FIG. 4. Top panel: Number of 4 He atoms as a function of 
the radial distance for the droplets of Fig. 3. Bottom panel: 
Number of 3 He atoms as a function of the radial distance for 
the same droplets. 



FIG. 1. Density profiles of 4 He728 + 3 Hejv 3 droplets for 
N 3 values from 1000 to 10000 in AA 3 = 1000 steps. To ease 
the figure, only the 4 He densities corresponding to a few A3 
cases have been plotted. Also shown is the density profile of a 
4 He728 drop immersed into liquid 3 He (dotted lines). Insert: 
bulk 3 He concentrations. The connecting solid line is to guide 
the eye. Also shown is the value corresponding to 4 He728 in 
liquid 3 He (dotted line). 



FIG. 2. Density profiles of 4 He 72 s + 3 HeiV3+ SF 6 droplets 
for for A3 values from 4000 to 10000 in A A3 = 1000 steps. 



FIG. 3. Bottom panel: Density profiles of 4 Hejv 4 + 
3 Heiooo + SF 6 droplets for N 4 = 35, 60 and 100. Top panel: 
Density profiles of 4 Hcjv 4 + 3 Heiooo + Xe droplets for the 
same A4 values. 
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